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MODIFIED  METHOD  OF  GEOMETRIC  ELECTROMAGNETICS 
FOR  THE  ANALYSIS  OF  RADIO  FIELD  IN  MARINE 
TROPOSPHERIC  WAVEGUIDES 

Anatoly  Bychkov,  Alexei  Bychkov 

Dept.  Radio  Engineering,  Sebastopol  Navy  Institute,  Sebastopol,  Ukraine 

The  interest  for  practical  use  of  tropospheric  waveguides  (TWG)  above  the  sea 
surface  is  still  high.  This  is  because  there  are  many  new  tools  for  remote  sensing  of  the 
atmosphere  boundary  layer  and  efficient  methods  for  calculation  of  radiowave  fields. 
Those  TWG  often  form  over  the  Black  Sea. 

An  efficient  method  for  calculation  of  radiofields  is  the  method  of  geometric 
electromagnetics  (GE)  based  on  the  Brillouin  conception.  It  has  been  used  to  determine 
radiofields  and  parameters  of  dielectrics  waveguides  and  microwave  devices  [1],  The 
method  GE  had  been  developed  for  waveguiding  systems  with  "rigid"  boundaries  of 
layered  environment.  However,  a sea  TWG  is  not  a system  with  "rigid"  boundaries.  The 
first  "wall"  of  TWG  is  the  surface  of  the  sea  with  parameters  £3,  |_t 3 , C3"  3 that  reflects 
incident  radiowaves.  The  layer  of  atmosphere  having  a lower  value  of  the  relative 
dielectric  permittivity  £2  will  be  the  second  "wall".  At  the  same  time  in  the  TWG.  up  to 

the  height,  h\,  the  permittivity  £|  varies  to  c2  smoothly.  If  a superposition  of  plane 

waves  (e.g.,  incident  and  reflected  waves)  satisfies  the  boundary  conditions,  then 
guided-wave  propagation  is  possible.  The  equation  for  the  resulting  field  in  TWG 
accounts  for  this  circumstance. 

The  Poynting  vector  77  of  a plane  radio  wave  propagating  in  TWG  (Fig.l)  forms 
angles  9 x , 0 y , 0 _ with  the  axes  of  coordinates: 

cos2Qx  + cos2  9 y + cos2  9.  = 1 . (1) 


At  the  height  of  TWG  h\  the  dielectric  permittivity  £j  starts  changing  smoothly 
and  reaches  e2at  the  height  h2,  hence  to  characterize  the  plane  radio  wave  in  the 
environment  with  a variable  index  of  refraction  we  can  use  the  equation. 
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where  cosQx,  cosQ  y,  cosQ  z are  projections  of  the  unit  normal  to  wave’s  front,  and 

n = a/s"  is  the  index  of  refraction.  Equation  for  a smoothly  varying  environment  is 

W(x,y,z)=  nW0(x,y,z)-  const , (3) 

where  function  W{x,y,z ) is  the  eikonal: 
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For  TWG  under  consideration,  the  function  W(x,y,z ) can  be  simplified  because 
cosQx  —0,  The  function  W(y,z)  might  be  presented  as  minimum  of  the  linear 
integral  jV?(iS  )dS  taken  from  S — 0 to  S = h\ . The  minimum  is  reached  if  the  curve 
S coincides  with  the  ray  trajectory  in  troposphere: 
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The  equation  (5)  is  used  for  determining  the  phase  of  radio  waves  reflected  from 
the  “upper  wall”  of  TWG.  A superposition  of  the  incident  wave  and  the  reflected  from 
the  media  interface  yields  the  following  electric  field  amplitude: 
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where  ky j — k\  • cosQ  y , kA  —k\  • cosQ Cp^/2  is  the  phase  of  the  coefficient  of 

reflection  from  the  “upper  wall”  of  TWG. 

The  equation  (6)  is  valid  under  the  conditions  of  total  interior  reflection  of  radio 
waves.  As  far  as  “upper  wall”  of  TWG  is  “soft”  (no  “rigid”  boundary  between  layers 
with  S j and  S2 ),  the  current  phase  of  the  coefficient  of  reflection  depends  on  the  angles 
of  incidence  of  radio  waves  and  on  the  variation  of  the  relative  index  of  refraction: 


Mt)=(VsW-1)-106+~-  = Wt)-i]- i°6+-^,  (7) 

kE  Ee 

where  RE  is  the  radius  of  Earth,  y is  a current  height  of  TWG  for  the  waves  with 

horizontal  or  vertical  polarization,  and  the  proper  phase  of  the  coefficient  of  reflection. 

To  account  for  the  vector  structure  of  the  field  it  should  be  transformed  to  new 
system  of  coordinates,  OX’ , 0 Yr , and  0 Z\  This  transformation  is  characterized  with 
the  Eiler  angles.  If  a radio  wave  has  vertical  polarization  and  jLij  = )U7  = 1 , 

a 1 ~ a 2 = 0 , then  we  can  write  the  final  equation  for  g>h  /2 : 
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where  «(y)  is  calculated  after  (7).  The  coordinate  varies  from  0 through  the 
inflection  point  of  the  profile  M(y),  i.e.  y — then  to  height  hj  that  is  the 
coordinate  of  curve  M{y ) where  the  atmosphere  is  standard. 

The  value  of  (p^/2  obtained  from  (8)  can  be  used  in  (6)  to  calculate  the  field  both 


inside  TWG  and  outside  it  at  the  heights  from  h\  to  /z2 . If  dielectric  permittivity  of 

troposphere  is  complex-valued  (conductance  a ^ 0 ),  then  the  equations  (5)-(8)  derived 
in  this  work  can  be  used  to  calculate  the  energy  losses  in  marine  TWG. 

The  modified  GE  method  that  has  been  discussed  in  this  work  enables  a more 
complete  account  of  TWG  effects  during  the  waving.  The  evaluation  of  a frontier 
conditions  by  the  “low  wall”  TWG  above  the  Black  Sea  shows  that  radio  waves  with 
length  about  10  cm  experience  almost  specular  reflection.  The  absolute  value  of  the 
reflection  coefficient  is  in  the  range  of  0.87  to  0.95.  The  phase  of  reflection  is  n for  the 
vertical  and  horizontal  polarizations  for  “weak”  and  “strong”  waveguides  of  evaporation 
having  h\  =10  m and  near-surface  TWG  having  h\  up  to  500  m when  the  sea  waving 
reaches  2 to  3 units  [2]. 

Modification  of  the  GE  method  allows  estimation  of  important  parameters  of 
wave  propagation  in  the  TWG:  wavelength,  maximum  phase  and  group  velocities, 
critical  angle  9cr  limiting  the  range  of  wave  propagation.  The  radiation  patterns  of 


radio  devices  must  be  oriented  into  angular  sector  that  is  determined  by  9cr , therefore 
knowledge  of  0C/.  has  great  practical  value  for  proper  feeding  of  TWG.  Our  theoretical 
results  have  been  confirmed  in  radar  experiments. 
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